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ID diagonally disordered chain with Frenkel exciton and long range ex- 
ponential intersite interaction is considered. It is shown that some states of 
this disordered system are delocalised contrary to the popular statement that 
all states in ID disordered system are localised. 

I. INTRODUCTION AND MAIN RESULT 

It is well known that all wave functions of translationary symmetric systems are delocal- 
ized. One of the most interesting properties of the homogeneous disordered systems is the 
possibility of localised wave functions. 

The mathematical problems of the theory of disordered systems are very complicated 
and for this reason the theory of disordered systems is not so well developed as the theory 
of symmetric systems. Despite this fact some statements related to disordered systems 
are considered to be well established and reliable. The above mentioned occurrence of the 
localised band is one of them. The next example of statement of this kind is that all states 
in ID disordered system are localised [4]. Recently appeared the reports [3,2] about the 
delocalisation in ID systems with intersite interaction in the form: J n ^ m = J/\n — m\ u , 1 < 
v < 3/2. In this letter we consider ID diagonally disordered chain with exponential intersite 
interaction and present arguments (computer simulations and theoretical treatment) in favor 
of partial delocalisation in this system. In this section we describe the system and present 
numerical results and in the next section we review the reasons which made us to study this 
system and present the approximate expression for the mobility edge. 

Let us consider ID Frenkel exciton in diagonally disordered chain. The mathematical 
problem is redused to the following random matrix of the Hamiltonian: 

H r y = e r 5 r y + w(r — r') r, r' — 1, N, (1) 

where 

w(r) — Vo exp — \r/R\, (2) 
Random values e r are supposed to be independent and having the distribution function: 

p(e) = | £ ^ ^ (Anderson's model) (3) 
I other cases 
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The thermodynamic limit JV — > oo is implied. To separate the localised and delocalised 
wave functions one should use some criterion of localisation. We use the number of sites 
covered by the wave function [7] determined as follows. Let us consider some eigen vector 
^ of the Hamiltonian (1) with components ^ r ,r = 1, JV. What contribution one should 
ascribe to the arbitrary site r? It is naturally to accept that this contribution is zero if 
\^ r \ 2 = and equal to unit if \^ r \ 2 =max{|\I/i| 2 , |\I/Ar| 2 }. So we come to the conclusion 
that the contribution of the arbitrary site r is |\I / r | 2 /| , I / | 2 na . r . The total number of sites n(\I / ) 
covered by normalised eigen function \& is the sum of contributions of all sites: 



n(¥) = £ 
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Delocalisation in (1),(2),(3) appear when R » 1. Below we study the properties of the 
eigen vectors of (1) with R = 20, v = 0.5, A = 4, JV = 1000. 

The dependance of number of sites covered by the wave function against corresponding 
energy for the Hamiltonian (1) is presented on fig.l (top). It is seen that n(E) is drastically 
increasing for energies higher than Eq ~ A. Additional calculations shows that n(E) do not 
depend on the number of sites JV for E < E and is ~ JV for E > E . For all these reasons 
we conclude that states below E are localised and states above E are delocalised. 

II. QUALITATIVE TREATMENT 

On our opinion the main properties of the above model which are responsible for the 
delocalisation are long range of intersite interaction R and the fact that function p{e) differs 
from zero only in the finite region. For these reasons for the qualitative interpretation we 
apply the following exactly solvable simple model of disordered system. Let the radius of 
interaction be infinite and write down the simplified Hamiltonian in the form: 

v 

H ry = 5 ry s r + — , r, r' — 1, JV (5) 

Taking advantage of the coherent potential approximation [5,4] one can show that the 
density of states for Hamiltonians (1) and (5) is coincide in the limit R — > oo, v — > 0, 2Rv = 
v. We show below that the Hamiltonian (5) has one delocalised and JV — 1 localised eigen 
functions. Consequently at least one delocalised function should appear in the set of eigen 
functions of the Hamiltonian (1) in the limit R — > oo. The desire to see how this take place 
was the starting point for our study of the Hamiltonian (1) with R » 1. Now let us turn 
to the proof of the above properties of the Hamiltonian (5). 

The equation for eigen vector e and eigen value A of the Hamiltonian (5) can be written 
in the form: 
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5 = Ee r (6) 



(6) gives an explicit expression for the eigen vectors of (5) as a functions of r and eigen 
number A. By substituting e r in the formula for S one can obtain the equation for the eigen 
values A: 

1^1 1 

^g— -f(A) = - (7, 

For Anderson's model (3) all quantities e r are differs from each other. For the graphical 
treatment of (7) the qualitative form of T(E) function is presented on fig.2. From fig. 2 one 
can see that N — 1 eigen values are belong to [0, A] and the last eigen value E m is not belong 
to this interval and in the thermodynamic limit can be determined from the equation: 



v ; J E-e A E-A v 

whence 



E m = z (9) 

1 — exp —(A/v) 

So in the thermodynamic limit E m is separated from any of e r by finite interval. From 
(6) one can see that sharp extremums of the wave function related to the localisation can 
appear if A G [0, A]. E m do not belong to this interval and we come to the conclusion that 
the corresponding eigen vector in the case of homogineous disorder is delocalised. Now let 
us show that all others eigen vectors are localised. For this reason introduce the Green's 
function in t-representation exp(?tH) r r which describe the dynamics of the wave function 
on the site r if it was equal to 1 on this site at t — 0. If the finite part of eigen states of the 
Hamiltonian H is delocalised this function goes down to zero when t — > oo and N — > oo. If 
Green's function do not decrease it means that the main part of eigen vectors is localised 
and the part of delocalised states is extremely small [4]. It is convenient to introduce the 
Green's function in E-representation: 

expitH= f G(E - i6) exp(iEt)dE 5^+0, t>0 (10) 

In the case of Hamiltonian (5) the Dyson's series for G: 

11111 11 

G rr , = b rr < + — v + — v 2 T(E) + (11) 

E — e r N E — e r E — e' r N 1 1 E - e r E - e' r 1 ; 

+i« 3 r(£) 2 - — + ... 

N 1 ' E-e r E-e' r 

can be exactly summed and give the following expression for G rr : 
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where 
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A E-A 



(13) 



In the thermodynamic limit the term ~ 1/N should be omitted and we come to the 
conclusion that the Green's function have a single pole E = e r . This corresponds to the 
oscillations of the wave function with constant amplitude and we can conclude that the main 
part of states are localised. It is easy to see that above described oscillations corresponds to 
the wave function localised on the site r and having an eigen value e r + 0(1/N). From fig. 2 
one can see that there are N — 1 eigen values of this kind and we come to the conclusion 
that the Hamiltonian (5) have iV — 1 localised states and one delocalised with eigen number 



Note that the appearance of the separated delocalised state for (5) is possible only if 
the distribution function p(e) is differ from zero in finite interval. For this reason we expect 
that delocalisation in (1) is also possible if p(e) is differ from zero in finite interval or at 
least goes down to zero rapidly enough. This statement confirms by calculations for Lloyd's 
model with p(e) = (l/7r)A/(A 2 + e 2 ) when no delocalisation was found. 

The energy dependance of number of covered sites for Hamiltonian (5) is presented on 
fig.l (bottom) for v = 20. Other parameters are the same as for the top picture. One can 
see that finiteness of the interaction radius results in appearance of delocalised states in the 
gap [A, E m ] but the boundary energy of spectrum and the mobility edge are the same for 
both Hamiltonians (1) and (5) and are equal E rn (9) and A respectively. 
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FIG. 1. Energy dependance of number of sites covered by the wave function for the Hamiltonian 
(1) with N = 1000, A = 4, w(r) = v exp-\r/R\, v = 0.5 , R = 20 (top). The same for the 
Hamiltonian (5) with v = 2Rvq = 20 (bottom). 
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FIG. 2. 



